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ABSTRACT f(t) 

The U.S. Navy and NASA are currently involved J 
in the design and development of an unsymmetric- 
skew-wing aircraft capable of 65° wing sweep and K 

flight at Mach 1.6. A generic skew-wing aircraft 
model was developed for 45° wing skew at a flight 
condition of Mach 0.70 and 3048 m altitude. At L, Li 

this flight condition the aircraft has a wing 
flutter mode. An active implementable control law 
was developed using the linear quadratic Gaussian M 
design technique. A method of modal residualiza- 
tion was used to reduce the order of the control- Pj 
ler used for flutter suppression. 
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SYMBOLS Q 

plant matrix Q ( s ) 

coefficient matrices of unsteady 
aerodynamic force approximation Qij(s) 

reduced-order controller state- 
space matrices R 


control and noise distribution 
matrices, respectively 

modal system input matrices 
corresponding to the low- and 
high-frequency parts of A 

lag coefficients 

state-space output matrix 

modal system output matrices 
corresponding to the low- and 
high-frequency parts of A 

generalized damping matrix 

forcing function in the s plane 


forcing function in time domain 

quadratic cost criterion 

regulator gain matrix, or general 
ized stiffness matrix 

robust and ordinary Kalman estima 
tor gain matrices, respectively 

generalized mass matrix 

pressure distribution of the jth 
mode 

state weighting matrix 

matrix of approximated aerodynami 
force coefficients 

the (i,j) element of a matrix of 
functions due to a mode 

vector of generalized coordinates 

control weighting matrix 

Laplace operator 

modal matrix 

input vector 

controller output 

gust velocity 

state vector and estimated state 
vector, respectively 

output vector 

mode shape of the ith mode 

control surface deflection and 
surface command, respectively 

random white noise excitation 



gust calibration factor 


a 

£ minimum singular values of the 

return difference matrix 

A similarity transform matrix 

Ai, A2 low- and high-frequency parts of 

A, respectively 

&)g wing gust input 

“m zero-mean white noise errors in 

the measurements 

INTRODUCTION 

Interest in oblique-wing aircraft designs has 
surfaced periodically since the 1940s. However, 
not until recently has the interest, technology, 
and mission of an oblique-wing design evolved into 
a full-scale flight research program. The U.S. 
Navy and NASA are currently in the design and 
development stage of implementing an oblique wing 
on an F-8 fuselage (Fig. 1) and evaluating the 
configuration to a maximum sweep of 65° and to 
Mach 1.6. 

The unsymmetric configuration and forward 
sweep of one semispan result in aeroelastic 
behavior distinctly different than that of 
straight, swept-back, or swept-forward wings. 

It should be noted that in addition to unsym- 
metric modeling characteristics, unsymmetric 
configurations will typically have signifi- 
cantly larger plant formulations since all 
degrees of freedom must be adequately repre- 
sented. Separation of an unsymmetric model 
into two smaller models (as is possible for 
symmetric and antisymmetric modes of a sym- 
metric aircraft) is not possible, because the 
response motion is coupled and not separable. 

To evaluate the analytical tools required 
for the analysis of an oblique-wing configuration, 
a generic skewed-wing model was developed. This 
model was used for the control system synthesis 
procedure described in this paper. The config- 
uration selected has a wing skew of 45° at a 
flight condition of Mach 0.70 and 3048 m altitude 
(a dynamic pressure of 23,892 N/m 2 ). 

Thi.s paper demonstrates the control synthesis 
design process required to develop a practical 
control law for stabilization of the flutter mode. 
This process involves 

1. formulation of the state-space model 
including independent wing actuators, a Oryden 
gust model [1], and s-plane approximations of 
unsteady aerodynamics; 

2. optimal full -state control law 
determination; 


3. robust output feedback control law 
determination; 

4. reduced-order (practical) control law 
formulation; and 

5. evaluation of the practical control law. 

MATHEMATICAL MODELS 

Aircraft Model 

The generic oblique-wing aircraft model used 
in the system synthesis process consists of a 
simple beam representation of the fuselage and 
wing. The structural model and the aerodynamic 
paneling required for the unsteady aerodynamics 
are represented in Fig. 2. 

The aircraft modal characteristics were devel- 
oped using NASTRAN analysis. At the selected 
sweep configuration (45°) and flight conditions 
(Mach 0.70, 3048 m altitude), the unaugmented 
aircraft has a flutter mode characterized as pri- 
marily wing bending but with some torsion. The 
m-vacuum mode shape characteristics of the wing 
alone are presented in Fig. 3 for the mode that 
is driven to the flutter condition with increasing 
dynamic pressure. 

Because the intent of this paper is to pre- 
sent a design synthesis process, the model order 
was reduced considerably; the final model con- 
tained a rigid-body (primarily pitch) mode along 
with three elastic modes. The model reduction 
process did not significantly affect the flutter 
mode characteristics. 

The formulation of the complete, integrated 
(structures, aerodynamics, and controls) state- 
space model for use in the analysis and design 
process follows that of Peele and Adams [2]. The 
aeroelastic equations of motion for the flexible 
aircraft can then be represented as 

Mcf + Dq + Kq = f(t) (1) 

where the matrices M, D, and K are the generalized 
mass, damping, and stiffness matrices, respec- 
tively, q is the vector of generalized coordi- 
nates, and f(t) is the vector of unsteady aero- 
dynamic forces; the dots denote differentiation. 
Transforming to the Laplace domain yields 

(Ms 2 + Ds + K)q(s) = F(s) (2) 

where s is the Laplace operator and F(s) is the 
aerodynamic forcing function in the s plane. The 
unsteady aerodynamic forces can be expressed as 

F(s) = Q(s)q(s) (3) 
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where the matrix Q(s) contains the generalized 
aerodynamic force coefficients. The individual 
elements of Q(s) are functions of both aircraft 
mode shapes and pressure changes resulting from 
motion in the various modes. The elements are 
defined as 

Qij(s) = jfjTz 1 (x,y)p J (x,y,s) dS (4) 

where z^x.y) is the mode shape of the 1th mode 
and Pj(x,y,s) is the pressure distribution of the 
jth reduced frequency determined from a lifting 
surface theory. Under the subsonic conditions 
relevant to this study, the unsteady aerodynamic 
force coefficients (elements of Q) were computed 
using the doublet lattice routine contained in the 
ISAC program [2,3]. Eight reduced frequencies 
were used, covering the range of 0 to 1.2 rad/sec. 
However, this procedure yields aerodynamic forces 
only for pure harmonic motion, and therefore only 
a finite number of frequencies can be selected. 

As a result, tabulated aerodynamic forces are 
expressed as a function of frequency. To apply 
modern control techniques, the tabulated aero- 
dynamics must be expressed in state-space form. 

If analytic continuity is assumed, the aerodynamic 
data can be expressed as a rational function 
approximation [4], such as. 


It is assumed that the actuators have suffi- 
cient power throughout the frequency range of 
interest and that aerodynamic hinge moments and 
inertial cross coupling do not affect control 
surface position. 

Gust Model 


The following second-order Dryden gust 
model [1] was also incorporated in the mathemat- 
ical model: 


0-273(1 + 4.114s) 
n (s + 0.421)2 


(7) 


where n is random excitation, a is used to cali- 
brate the gust intensity (to 1 ft/sec in the cur- 
rent example), and vg is the output gust velocity. 

State-Space Equations 

The design model is obtained by combining the 
aircraft (including the linearized form of the 
unsteady aerodynamics), the actuators, and the 
gust model dynamics and can be represented in the 
state-space equation form as 


x = Ax + B^u + B 2 ug 


( 8 ) 


n 

Q(s ) = Ao + Ais + A 2 S 2 + X A 4+ 2 [s/(s + b A )] 

t=l 

(5) 

where Ag to A n+ 2 are coefficient matrices of 
unsteady aerodynamic force and b* are "lag" coef- 
ficients. A least squares approach can then be 
used to determine the matrices Ao, Aj., A 2 , .... 
A n +2* The lag coefficients bj, b 2 , .... b n are 
selected specifically for the analysis; the number 
of lag terms and their values are important in 
obtaining good approximations of the tabulated 
aerodynamics. In the oblique-wing design model 
described in this paper, two lag terms were used 
in generating the s-plane fit of the unsteady 
aerodynamics. A typical fit of the approximation 
to the tabulated data for one element of the Q(s) 
matrix is presented in Fig. 4. 

Actuator Model 

Left and right wing actuators were modeled 
independently because the synthesis process deter- 
mines unique control laws for each surface. The 
following third-order model relates the control 
surface deflection 6 a to the control surface com- 
mand 6 3c 

5 a 54,080 (6) 

5 a c (s+20)[s 2 + 2(0.7)52s + 522] 


y = cx + ofo (9) 

where x is the state vector, ug the wing gust 
input (unit white noise), aim the measurement 
noise, u the control input vector (2 x 1), and y 
the measurement vector (3 x 1); A, Bi, B 2 , and C 
are plant equation, control, noise distribution, 
and state-space output matrices, respectively, of 
suitable dimensions. The state vector contains 
24 states, including the rigid body mode, flexible 
mode deflections, flexible mode rates, unsteady 
aerodynamic states, actuator deflection and rate 
states, and wind gust states. Eight states result 
from the structural modes retained, eight from the 
two-lag-term set of approximated unsteady aerody- 
namics, six from the two actuators, and two from 
the gust model. The three outputs are the acce- 
lerations at the center of gravity, right wingtip, 
and left wingtip. 

CONTROL LAW DESIGN 
Optimal Controller 

The linear quadratic Gaussian (LQG) method 
is viable for the design of multi-input multi- 
output controllers. The active control syn- 
thesis is based on LQG theory but is modified 
to accommodate the high-order model of the 
aircraft [5,6]. The design process involves the 
following steps state-space model generation; 
full-state feedback design, estimation of states 
from available measurements, and development of 
reduced-order controller. 
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The state-space model of the aircraft is 
defined by Eqs. (8) and (9). A full-state feed- 
back control law. 


z = A z + T _1 Ly 
u c = -KTz 


u = -Kx 


( 10 ) 


where T is the modal matrix and 


(15) 

(16) 


is determined by minimizing a quadratic cost func- 
tion [7] 


J 



u T Ru) 


dt 


( 11 ) 


where Q and R are suitable weighting matrices. 

Because direct measurement of all states of an 
aeroelastic system is not feasible, it is necessary 
to estimate states from available measurements. A 
Kalman filter is used for estimation of the states. 
The estimator dynamics are given by 


i = (A - BiK - LiC)x + Liy (12) 


where l_i is the Kalman estimator gain matrix and x 
is the vector of estimated states. However, 
systems designed using a Kalman estimator are con- 
ditionally stable, have poor gain and phase mar- 
gins, and have high bandwidth [8]. The input 
noise procedure of Doyle and Stein [9] can be used 
to synthesize a robust Kalman estimator. This 
procedure involves compromising root-mean-square 
(rms) response activity against robustness. The 
optimal controller, which is of the same order as 
the aircraft model used for synthesis, and the 
controller output u c are defined as 

‘x = (A - BK - LC)x + Ly (13) 


A = T - 1 A au gT 


A au g = (A - BjK - LC) 
x = Tz 

This can be expanded to 


' • m 
Z 1 


Al 

0 ’ 

Z 1 

+ 

V 

z 2. 


0 

A 2. 

_ z 2. 


_ b 21_ 


u c = [Cl C2] 


Z 1 

L Z 2J 


(17) 


(18) 


where Aj, Ag, Bij, B 21, Cj, and C2 are the matri- 
ces corresponding to the low- and high-frequency 
parts of the original system matrices. In the 
low-frequency portion the dynamics are retained, 
while in the high-frequency portion only the sta- 
tic terms with zero response time assumed for the 

dynamics are retained. Setting £2 = 0, 


z 2 = '^2 B 21 y 

Substituting Eq. (19) into Eq. (18), 


In 


= C^Zj - C 2 (A 2 ® 2 jy) 


(19) 


( 20 ) 


u c = -Kx 


(14) 


Rewriting the reduced system. 


where L is the robust Kalman estimator gam matrix. 


z x = A^ + 8 n y = Azj + By 


( 21 ) 


Practical Controller 


The full-order optimal controller consisting 
of a robust Kalman estimator together with opti- 
mal state feedback gains imposes an unnecessarily 
large implementation cost. A reduced-order con- 
troller that approximates the full-order optimal 
controller can be found that imposes little deg- 
radation in performance [10]. A modal residual - 
ization technique [11] can be used to reduce the 
order of the controller. An attempt is made to 
approximate the full-order controller with a 
lower-order approximation while maintaining the 
desired characteristics of the original control- 
ler. A similarity transform A is employed on the 
full -order controller described by Eqs. (13) and 
(14) to obtain 


-In 


u c ~ 6^Zj _ CjjAg ®21^ = Cz - Dy 


where 


A = A, 


B = B 


11 


C = C, 


D = -C 2 A 2 B 21 


( 22 ) 
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The eigenvalues of the reduced system are the 
eigenvalues retained in the Aj portion. A block 
diagram of the plant and the reduced-order 
controller is presented in Fig. 5. 

FLUTTER SUPPRESSION APPLICATION 

The control law synthesized by the method 
outlined in the previous section was applied to 
the design of an active flutter suppression 
controller for an oblique-wing aircraft. A 
generic 45°-wing-skew structural model was devel- 
oped to simulate flutter at a subsonic flight 
condition of Mach 0.70 and an altitude of 3048 m. 
The unstable eigenvalue pair at this flight con- 
dition (0.50 ± j 14 . 37 ) represents primarily wing 
bending with some torsion. 

The design objective was to stabilize the 
aircraft without exceeding the specified rms 
control activity so that saturation would not 
occur. Based on actuator limitations, the rms 
deflection of the aileron was limited to 5° and 
the deflection rate to 30 deg/sec. In addition 
to stabilizing the aircraft with low surface 
activity, it is required that the controller be 
robust. The controller considered here is 
multi-input, multioutput: The right and left wing 

control surfaces are independent of each other 
because of the unsymmetric nature of the aircraft. 
Robustness of the multiloop control system is 
evaluated by using the singular values of the 
return difference matrix [12,13,14]. 

Linear quadratic Controller Design 

A full -order control law, identical to the LOG 
solution, was obtained first for comparison pur- 
poses. Initial values of the weighting matrices Q 
and R of Eq. (11) were selected as null and iden- 
tity, respectively [6], All stable eigenvalues 
remain unchanged and all unstable eigenvalues are 
rotated about the imaginary axis [7]. With the 
resulting values of the full-state feedback gains, 
the remaining design process was executed (robust 
output estimation and controller reduction). 

Though the rms control activity of the reduced 
controller was within specified limits, the con- 
troller was not very robust. To improve the 
robustness of the final reduced controller, while 
still retaining low surface activity, parametric 
variations of Q and R were performed. A matrix R 
with values of 50,000 along the diagonal and a 
positive-definite Q with values of 0.00001 along 
the diagonal gave satisfactory characteristics. 

The minimum singular value of the return differ- 
ence matrix, jj, for full -state feedback was always 
>1, as shown in Fig. 6, the rms control values are 
given in Table 1. 

Full-Order Controller 

Because all the states are not available for 
feedback, it is required that all states be esti- 
mated so that regulator gains can be used. The 
controller size is the same as that of the 
aircraft plant, and will be referred to as the 
full-order controller. 


A Kalman estimator is used to estimate the 
states. The rms values for the control activity 
are given in Table 1. The robustness, as indi- 
cated by the minimum singular value plot of Fig. 7, 
is relatively poor, which is characteristic when 
using an ordinary Kalman estimator [8]. A design 
procedure described by Doyle and Stein [9] is used 
to improve the robustness of the Kalman estimator. 
This method involves applying extra process noise 
to the control input of the aircraft during esti- 
mator design. Figure 8 shows the minimum singular 
value plot for the return difference matrix for 
the estimator designed using the additional noise 
(robust Kalman estimator). The improved stability 
margin does increase the rms control activity (as 
shown in Table 1), but it is still within the 
specified limits. 

Reduced-Order Controller 

The robust Kalman estimator, together with 
the optimal feedback gains, constitutes an opti- 
mal implementable controller. It is, however, 
impractical to implement this controller because 
of the cost involved. The cost of implementation 
can be reduced by developing low-order approxima- 
tions to the optimal controller; approximations 
are referred to as practical or reduced-order 
controllers. The reduced-order controller must 
achieve closed-loop stability, have satisfactory 
control activity, and be robust. 

A seventh-order controller was obtained by 
using the method of modal residual ization. Table 2 
shows the eigenvalues of the full-order robust 
controller and the eigenvalues retained in the 
reduced-order controller. Figure 9 shows the step 
response of the full -order and reduced-order 
controllers and the extent of degradation in the 
response due to reduction. 

The rms control activity for the reduced- 
order controller is shown in Table 1, and 
Fig. 10 is a plot of minimum singular values 
for this case. Even though there is some 
degradation caused by controller order reduc- 
tion, the rms control activity and the stability 
margin are considered acceptable. 

CONCLUSIONS 

An implementable flutter controller for a 
45°-skew oblique-wing aircraft mathematical model 
was designed using the LQG design methodology. 
Kalman estimators produced low stability margins, 
however, the Doyle-Stein procedure for robust 
estimator design can be used to improve these 
margins to acceptable values without excessive 
surface activity. A modal residualization tech- 
nique was used to obtain a reduced-order con- 
troller that satisfied the performance require- 
ments and can be implemented. 

New controllers will be designed as improved 
models (based on the actual flight configuration) 
of the oblique wing are made available. Actual 
implementation may require gain scheduling as a 
function of wing skew as well as other parameters. 
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Table 1 RMS responses at flutter conditions 



Right 

wing 

Left 

wing 


• 

6, deg S 

, deg/sec 

5, deg 

6, deg/sec 

Full -state feedback 

1.58 

7.91 

0.28 

3.58 

Full-order controller 
with Kalman estimator 

1.58 

12.58 

0.41 

5.64 

Full-order controller 
with robust Kalman 
estimator 

2.06 

11.21 

0.41 

5.02 

Reduced-order 

controller 

1.51 

10.04 

0.41 

4.90 
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Table 2 Full -order and reduced-order 
controller eigenvalues 


Full-order 
control ler 
eigenvalues 

Reduced-order 
controller 
ei genvalues 

-0.2531 

+ 

0.00001 

-0.2531 + 

O.OOOOi 

-0.4578 

+ 

0.0845i 

-0.4578 + 

0.0845i 

-0.4578 

- 

0.08451 

-0.4578 - 

0.08451 

-6.3566 

+ 

0.00001 



-2.5815 

+ 

6.36171 

-2.5815 + 

6.3617i 

-2.5815 

- 

6.36171 

-2.5815 - 

6.3617i 

-0.3773 

+ 

13.1104! 

-0.3773 + 

13.11041 

-0.3773 

- 

13.11041 

-0.3773 - 

13.11041 

-3.4330 

+ 

15.36501 



-3.4330 

- 

15.36501 



-20.0399 

+ 

O.OOOOi 



-28.5469 

+ 

9.70161 



-28.5469 

- 

9.70161 



-36.0773 

+ 

0.3261i 



-36.0773 

- 

0.32611 



-37.0844 

+ 

2.44341 



-37.0844 

- 

2.44341 



-40.1968 

+ 

1.37961 



-40.1968 

- 

1.37961 



-36.1567 

+ 

37.22251 



-36.1567 

- 

37.2225i 



-53.6588 

+ 

O.OOOOi 



-34.5946 

+ 

44.06571 



-34.5946 

- 

44.06571 





Fig. 1 Oblique-uing configuration. 
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Fig. 2 Generic model (aero panels and node points). Vertical 
fin shown in X-Y plane. 


Wing 

displacement 



Wing 

rotation 



wingtip wingtip 


Fig. 3 In-vacuum flutter mode shape character- 
istics at 2.239 Hz. 45 0 skew. 


Imaginary 



Fig. 4 Unsteady aerodynamics comparison. 



B 2% D0 > m 



Fig. 5 Cloaed-loop system with reduced-order 
controller. 




Fig. 6 Minimum singular values for full-state Fig. 7 Minimum singular values for full-order 

feedback. controller with Kalman estimator. 


Minimum 
singular value, 
a 




Time, sec 


Fig. 8 Minimum singular values for full-order 
controller with robust Kalman estimator. 


Fig. 9 Step response comparison of full-order 
controller and reduced-order controller. 


Minimum 
singular value, 

o 



Fig. 10 Minimum singular values for reduced- 
order controller. 
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